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U q (l ) = £(q) (q> 1) , U q (t) = t + 0(t 2 ) , (15.46) 

Notar que Li q {t) diverge en t = 1 para q < 1. 

Volviendo al cálculo de /i, vemos que la (15.43) se convierte en 
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n = — / -+ vi 2mnkT \ U 3/2 (e^) 

\ h 2 ) 3/2 


(15.47) 


la cual, recordando que 


12jimkT c ^ 3/2 n 
\ h 2 ) ~ V£(3/2) 


(15.48) 


se escribe finalmente como 


IT \ 3 ' 2 LÍ3, 2 ( e ft‘) _1_ 

T c ) ^(3/2) ’ 1 e-^-1 


(15.49) 


Para un sistema macroscópico n es un número enormemente grande (del orden del número de 
Avogadro), de modo que n { /n es despreciable frente a la unidad, para todo T >T C ; por lo tanto, 
tenemos que, con buena aproximación 


j (T\ 3 , 2 U 3 / 2 (e^) 

UcJ e(3/2) 


(T > T c ) 


(15.50) 


Invirtiendo esta relación obtenemos finalmente el potencial químico (Fig. 15.2). Se puede ver 
(no damos aquí los detalles) que en el intervalo 0 < T < T c , es negativo y no se anula nunca 
(salvo para T = 0); su valor absoluto es muy pequeño, siendo del orden de kT / n. 



Fig. 15.2. Potencial químico de un gas de Bosones no interactuantes. Con la línea de trazos 
se muestra el valor clásico de ¡u, dado por la ec. (15.51). 
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El límite e^ 1 « 1 de la (15.50) permite obtener la expresión de [i para un gas perfecto clásico. 
Usando la segunda de las (15.46) obtenemos 

Los números medios de ocupación 

Es interesante calcular los números medios de ocupación, dados por la (14.28). La Fig. 15.3 
muestra ñ(e) para diferentes temperaturas del gas. El gráfico es semilogarítmico pues así se 
aprecian mejor los apartamientos desde la distribución clásica de Maxwell-Boltzmann, que se 
representa por medio de rectas de pendiente 1 / kT. 


úciásico kT\ n. 


1 


£(3/2)1 T c 



e/kT c 

Fig. 15.3. Número de ocupación medio en función de la energía del estado para un gas de 
Bosones. Las curvas corresponden a T/T c = 1,1.5, 2, 3, 4 y 5. 

Se puede observar que las curvas para TI T c = 3, 4 y 5 son prácticamente rectas, lo que indica 
que el comportamiento del gas para esas temperaturas es clásico, como es de esperar porque en 
esos casos ñ(e) es siempre mucho menor que la unidad. En cambio, las curvas correspondientes 
a TIT C = 1,1.5 y 2 se apartan de las rectas clásicas para s por debajo de 2 - 2.5 kT c , donde ñ(e) 
ya no se puede despreciar frente a la unidad. 


El calor específico de un gas de Bosones 

Vamos a estudiar ahora el calor específico de un gas de Bosones. Recordando que nk = NR , 
donde N es el número de moles, tenemos que el calor específico molar está dado por 


~ _ R_(dE\ 

Cv ~~ñk\fff) v 


(15.52) 
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Por lo tanto para determinar c v tenemos que calcular primero la energía interna del gas, dada 
por 


E = y e¡n ■ = y „ £ ' . -s í 


ef(e)de 


= 2 jtV 


(2m\ 312 í e 3,2 de 


e P(e-n) ^ Jlr \h2) I e /3(e-,w) _ \ 

0 0 


(15.53) 


Haciendo el cambio de variable fíe = z y usando la (15.35) y la (15.44) resulta 

E = - nkT [—\ U 5l2 (e^) 

2£(3/2)\T c ) 5/2 


(15.54) 


Aquí conviene considerar por separado los casos T <T C y T >T C . 

Cuando T <T C , podemos suponer que f/»0y por lo tanto Li 5/2 (e^ u ) = Li 5/2 (1) = £(5/2). Re¬ 
sulta entonces que 


E = nk- ^ 5I2) T l' 2 (/ < /;) 

2 £(3/2 )T c 3/2 c 


(15.55) 


y por lo tanto 


3/2 


Cy — —R 


, (t < t c ) 
4 £(3/2 )[T C J 


(15.56) 


Cuando T >T C , obtenemos de la (15.52) 


c v = 


3 R 


(L 


3/2 


2£(3/2){T c J 2 5/ - 3/2 


E 1 Id [i 
kT~k\dT 


(15.57) 


Para obtener la expresión de (d[i/dT) tl derivamos la (15.50). Resulta 


kT 


_UdE\ 3 Li 3/2 (e^) 
k\¿rr) n 2U m (e^) 


(15.58) 


Obsérvese que tanto ¡i como (dfj,/ dT) n tienden a cero para T -» T c pues Li 1/2 (í) diverge en 
t = 1. Usando la (15.58) obtenemos finalmente 


3/2 


C V 


3 * ^nAi 5/ ,(^)- 3Li ^ )2 i 


2£(3/2)\T c ) 2 


2 Li 


1/2 


(e^) J 


(r>r c ) (15.59) 


Cuando T = T C tanto la (15.56) como la (15.59) convergen a 

c v (T c ) = —R ^ (5/2) = 1.92567 R 
V c 4 £(3/2) 


(15.60) 


En la Fig. 15.4 se muestra el comportamiento de c v . 
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Fig. 15.4. Calor específico molar a volumen constante de un gas de Bosones. La línea de 
trazos indica el valor clásico 3R/2. Para T > T c el calor específico del gas de Bosones es 
mayor que el valor clásico, y aumenta al disminuir T hasta alcanzar en T c el máximo dado 
por la ec. (15.60). Al disminuir la temperatura por debajo de T c , el calor específico del gas 
degenerado disminuye rápidamente con T. En T = T c la derivada de y / cfT es discontinua. 

Es interesante comparar las propiedades que acabamos de estudiar con el comportamiento del 
4 He líquido a bajas temperaturas. Como el 4 He tiene spin cero cumple la estadística de Bose- 
Einstein y por ser un gas inerte, las tuerzas interatómicas son las de van der Waals, que son muy 
débiles. Por ese motivo el helio se licúa a una temperatura muy baja (el punto normal de ebulli¬ 
ción del 4 He es de 4.2 °K), y sólo se solidifica bajo presiones muy grandes. Al enfriar el 4 He lí¬ 
quido en contacto con su vapor, sus propiedades cambian bruscamente a T = 2.17°K. Para 
T > 2.17 °K se comporta como un líquido normal, llamado helio I. Para T < 2.17 °K presenta 
propiedades muy llamativas, por ejemplo fluye por capilares muy finos como si no tuviese vis¬ 
cosidad. Esta forma se llama helio II, y sus características se describen por medio del modelo de 
dos fluidos, que trata el helio II como una mezcla de un fluido normal y otro superfluido, sin in¬ 
teracción viscosa entre sí. El fluido normal tiene todas las propiedades familiares de los fluidos. 
El superfluido, en cambio, tiene características muy curiosas: su entropía es nula, y no tiene vis¬ 
cosidad. El hecho que en un gas de Bosones degenerado hay dos clases de partículas cuyo com¬ 
portamiento es muy diferente sugiere que las propiedades del helio II se relacionan con el fe¬ 
nómeno de la degeneración. Hay muchas analogías entre el helio II y el gas de Bosones degene¬ 
rado, además de las que vamos a comentar. En la transición helio I-helio II se observa un com¬ 
portamiento anormal del calor específico (Fig. 15.5). Como la curva experimental recuerda la 
letra griega A (lambda), el paso de helio I a helio II se denomina transición A, y la temperatura 
de transición se llama punto A. La semejanza entre las Fig. 15.4 y 15.5 llevó a London 9 a inter- 


9 Fritz Wolfgang London y Walter Heitler desarrollaron en 1927 el primer tratamiento cuántico de la molécula de 
hidrógeno. F. London y su hennano Fleinz formularon en 1935 le teoría fenomenológica de la superconductividad, 
un fenómeno estrechamente relacionado con la superfluidez. 
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pretar que la transición A del 4 He es una manifestación de la condensación de Bose-Einstein. 
Usando la ec. (15.35) calculó la temperatura de condensación de un gas perfecto de átomos de 
helio, para una densidad igual a la del helio I en el punto A experimental. Con V/n = 27.6 
cm 3 /mol obtuvo T c = 3.13 °K, lo cual no está demasiado lejos del valor medido T x = 2.17 °K. 



Fig. 15.5. Comportamiento del calor específico del 4 He, tal como resulta de las medicio¬ 
nes. Se puede apreciar la transición A a T = 2.17 °K. 

Hay diferencias importantes entre el comportamiento del helio II y el gas de Bosones. Entre ellas 
se cuenta que la dependencia en la presión de la temperatura de transición es distinta. Además, 
las mediciones más precisas muestran que en el punto A, c v no se mantiene acotado como ocu¬ 
rre para el gas de Bosones, sino que tiene un infinito logarítmico. Si embargo estas diferencias 
no deben sorprender, pues la teoría que hemos desarrollado corresponde a un gas perfecto, y no 
se puede aplicar tal cual al helio líquido, en el cual hay interacciones interatómicas, que (evi¬ 
dentemente) no se pueden ignorar. Por consiguiente es razonable concluir que la transición A del 
4 He es lo análogo para un líquido de la condensación de Bose-Einstein de un gas perfecto. Esta 
conclusión se refuerza si se tiene en cuenta que el 3 He, cuyos átomos tienen spin 1/2 y por lo 
tanto son Fermiones, no presenta transición 10 en el punto A del 4 He. 

La superconductividad (que no tratamos por razones de espacio) es un fenómeno análogo a la 
superfluidez. Según la teoría BCS (formulada en 1957 por John Bardeen, León N. Cooper y 
John R. Schrieffer), a temperaturas muy bajas los electrones de conducción se agrupan de a pa¬ 
res {pares de Cooper) formando una suerte de cuasimolécula. Como los pares se comportan 
como Bosones, se pueden condensar y esto explica las propiedades de los superconductores. 


10 Para el 3 He, la superfluidez se observa a temperaturas inferiores a 3x1 (T 3 °K. Esto se debe a que a esas 
temperaturas tan bajas los átomos del 3 He se unen formando pares, de igual modo que los electrones en un 
superconductor. Los pares tienen spin entero y por lo tanto obedecen a la estadística de Bose-Einstein. Tan pronto 
se forman los pares, éstos sufren una condensación, y por lo tanto todo lo dicho para el 4 He se aplica también al 3 He 
cuando éste se encuentra a temperaturas para las cuales se han formado los pares. 
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Durante muchos años se careció de evidencia directa de la existencia de la condensación de 
Bose-Einstein, aunque se creía que ocurría en el helio líquido. La situación cambió en 1995 
cuando Michael H. Anderson, Jason Enscher, Michael Matthews, Cari Wieman y Eric Cornell, 
observaron la condensación de Bose-Einstein en un gas de 87 Rb con una densidad de 2.5xlO ls 
átomos/m 3 , enfriado por debajo de 0.17x10 6 °K. El gas estaba confinado por un campo magné¬ 
tico, y para detectar la condensación se desconectaba dicho campo y poco después se registraba 
la densidad espacial de las partículas. Aquellas cuya velocidad es apreciable se alejan rápida¬ 
mente de la región de confinamiento, pero los átomos del condensado, cuya cantidad de movi¬ 
miento es nula, permanecen en dicha región. Si hay muchas partículas en ese estado, se esperaba 
ver un pico de la densidad en la región donde el gas había estado confinado. Esto fue efectiva¬ 
mente observado, pero sólo si la temperatura está por debajo de T c . Al enfriar el gas por debajo 
de T c se encuentra que el pico se refuerza, tal como está previsto por la teoría. Ese mismo año, 
en forma independiente, Wolfgang Ketterle y sus colaboradores realizaron experimentos simila¬ 
res pero algo más sofisticados con átomos de sodio y de litio, y también observaron la condensa¬ 
ción de Bose-Einstein. Por estos trabajos Cornell, Ketterle y Wieman fueron galardonados con el 
Premio Nobel para la Física en 2001. 

El gas de fotones y la radiación de cuerpo negro 

Vamos a estudiar ahora la radiación de cuerpo negro, considerándola como un gas perfecto de 
fotones en equilibrio térmico. De acuerdo con la hipótesis de Einstein (ver el Capítulo 4) vamos 
a suponer que los fotones son partículas de masa nula y energía dada por la ec. (4.16), esto es 

£¡ = hv¡ (15.61) 


Como los fotones tienen spin 1, son Bosones y corresponde aplicar la estadística de Bose- 
Einstein. Al usar la distribución (15.24) hay que tener en cuenta que los fotones pueden ser emi¬ 
tidos y absorbidos por las paredes de la cavidad (este es el mecanismo que garantiza el equilibrio 
de la radiación de cuerpo negro a Ty V fijos), y entonces n no está determinado a priori, a dife¬ 
rencia del gas de Bosones que tratamos antes. En este caso n está determinado por la condición 
de equilibrio a T y V fijos 11 . Esto implica que el potencial químico del gas de fotones es idénti¬ 
camente nido. En efecto, para todo cambio espontáneo T y V fijos se cumple 

AF < 0 (15.62) 

donde F es la función de Helmholtz. Luego en el equilibrio F es mínimo, lo que lleva a la condi¬ 
ción necesaria (dF / dn) TV = 0. De aquí se obtiene entonces 

p = (—) = 0 (15.63) 

\dn) TV 

Teniendo en cuenta este hecho la distribución (15.24) queda de la forma 

ñ { = —¡r^ - fotones (15.64) 

e P E i ~ i 


11 ver el Capítulo 8 de Termodinámica e Introducción a la Mecánica Estadística. 
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A partir de la (15.64) es muy sencillo deducir la distribución espectral de Planck. En efecto, re¬ 
cordemos que la densidad de estados por unidad de intervalo de frecuencia es (ec. (4.1)): 

A(v) = ^v 2 (15.65) 

c 3 

Por lo tanto la densidad de energía ü(v,T ) (energía por unidad de volumen en el intervalo de 
frecuencia entre vy v + clv) está dada por 

ü(v,T) = V~ l N{v)ñ{e)e (15.66) 


Sustituyendo (15.63) y (15.64) en (15.66) se obtiene la distribución espectral de Planck (ec. 
(4.12)). 


8 Jtv 2 hv 
u(v,T)~ hv/kT 


c J e 


-1 


Es fácil verificar que el máximo de u(v,T ) ocurre para 


kT 

v = v m = 2.82144 — 
h 


(15.66) 


(15.67) 


La (15.67) expresa la Ley de desplazamiento de Wien, que fue obtenida originalmente en la 
forma v m ~ T en base a un argumento puramente termodinámico. 

Integrando la (15.66) sobre todas las frecuencias obtenemos la densidad de energía radiante 


k-ju(v,T)dv e ,J tT _ í dv 

0 


8 Jl 5 k 4 
\5c 3 h 3 


(15.68) 


La (15.68) es la Ley de Stefan-Boltzmann p EtV = (4 o / c)T 4 , y podemos entonces expresar la 
constante de Stefan-Boltzmann como 


2jt 5 k 4 
1 5c 2 h 3 


(15.69) 


Es oportuno comentar brevemente el significado físico de la distribución (15.64) para aclarar 
algunos aspectos de la radiación mencionados en capítulos anteriores. La Fig. 15.6 muestra la 
distribución de Planck (15.66) y la Fig. 15.7 los números de ocupación medios (ec. (15.64)) por 
medio de curvas universales, obtenidas adimensionalizado le energía de los fotones con el factor 
kT y ü(v,T) con el factor aT 3 , donde a = 60 ho I it 4 kc. Para comparar, mostramos también en 
esas figuras la posición de v m , y los resultados de la teoría clásica de Rayleigh-Jeans (líneas de 
trazos) y de la fórmula empírica de Wien mencionada en el Capítulo 4 (líneas de puntos). 

Por de pronto, comparando las Figs. 15.3 y 15.7 vemos que el comportamiento de ñ es análogo 
al de un gas de Bosones a T = T C (pues en ambos casos /x = 0). Por lo tanto, cualquiera sea su 
temperatura, el gas de fotones se comporta de manera semejante a un gas de Bosones a T = T c . 


12 Ver el Capítulo 15 de Termodinámica e Introducción a la Mecánica Estadística. 
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hv/kT 

Fig. 15.6. La distribución de Planck. Representamos también la distribución de Rayleigh- 
Jeans (línea de trazos), que reproduce para bajas frecuencias el resultado de Planck, y la 
distribución de Wien (línea de puntos), que tiene el comportamiento correcto para frecuen¬ 
cias altas. 



Fig. 15.7. Número medio de ocupación de los estados de fotones para la radiación de 
cuerpo negro. Las líneas de trazos y puntos corresponden a las fórmulas de Rayleigh-Jeans 
y de Wien, respectivamente. 

Comparando las Figs. 15.6 y 15.7 notamos que el ámbito de validez de la fórmula de Wien co¬ 
rresponde al rango de frecuencias para el cual ñ « 1, esto es, el rango en el cual los fotones 
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cumplen la condición (15.11) y entonces se comportan como partículas clásicas (esto es, distin¬ 
guibles). En efecto, la fórmula de Wien no es otra cosa que la distribución clásica para un sis¬ 
tema de partículas de masa nula y energía e = hv . Este rango, que podríamos llamar entonces 
de partícula clásica, abarca las frecuencias mayores que « 2kT/ h y comprende un 80% de la 
energía radiante total. Para comprender lo que esto implica daremos algunos valores. La radia¬ 
ción electromagnética que proviene de casi todas las fuentes naturales es de origen térmico 13 y 
su distribución espectral es, grosso modo, del mismo tipo que la del cuerpo negro. La luz solar 
corresponde a 5700 °K, y todos los fotones de longitud de onda más corta que la del cercano in¬ 
frarrojo están en el rango v > 2 kT/ h; para una temperatura de 300°K, lo están todos los fotones 
de longitud de onda menor que unos 25 pm. Por lo tanto, del infrarrojo en adelante, los fotones 
provenientes de fuentes naturales se comportan como partículas clásicas, y se describen me¬ 
diante paquetes o trenes de ondas que no se solapan; se trata entonces de radiación incohe¬ 
rente 14 . En estas condiciones los fenómenos de interferencia y difracción a que dan lugar se de¬ 
ben al comportamiento de cada fotón individual (cuyo paquete de ondas, por ejemplo, al pasar 
por las rendijas del experimento de Young interfiere consigo mismo). Asimismo, estos fotones 
interactúan con la materia de a uno por vez, como vimos en el Capítulo 4 al tratar el efecto fo¬ 
toeléctrico, el efecto Compton y otros fenómenos. 

En cambio la Ley de Rayleigh-Jeans vale para frecuencias bajas, para las cuales ñ » 1. En este 
rango el conjunto de los fotones que ocupan cada estado de energía e¿ = hv i está descripto por 
un paquete de ondas que es una superposición coherente de los paquetes de onda individuales. 
Por este motivo, cuando ñ es muy grande los campos eléctricos y magnéticos del conjunto tie¬ 
nen las propiedades de una onda clásica, caracterizada por la fase, además de la frecuencia y la 
amplitud. Puesto que en este dominio los cuantos interactúan con la materia en conjuntos cohe¬ 
rentes, es decir colectivamente y no individualmente, la noción misma de fotón pierde significa¬ 
ción. De hecho, en este régimen los fotones son inobservables. El rango de onda clásica 15 abarca 
las frecuencias por debajo de ~ 0.01 kT/ h y comprende menos del 3% de la energía radiante to¬ 
tal. Para 5700 °K corresponde a A > 0.025 cm y para 300°K, A > 0.5 cm. 

Los dos comportamientos que se acaban de comentar representan los casos extremos de altas y 
bajas frecuencias. La radiación térmica de frecuencias intermedias (O.OlkT/h < v < 2kT/h) 
tiene parte de las características de onda clásica y parte de las características de partícula clásica, 
ya que la transición entre ambos regímenes es continua. Vemos así que la descripción cuántica 
de la radiación incorpora la dualidad onda-corpúsculo de la cual hablamos en el Capítulo 4, y le 
da un significado más preciso. 

Para terminar esta discusión de la radiación de cuerpo negro, observemos que si bien la presente 
deducción de la (15.66) y la que dimos en el Capítulo 4 son esencialmente equivalentes, hay una 
diferencia conceptual entre ellas que conviene subrayar. Mientras en el Capítulo 4 consideramos 


13 Hay algunas excepciones, como por ejemplo la radiación de sincrotrón emitida por fuentes cósmicas 
(radiogalaxias y quasars). 

14 Aquí estamos empleando los términos “coherente” e “incoherente” en un sentido algo diferente del habitual en la 
Óptica clásica. 

15 Puesto que éste es el rango donde se ponen de manifiesto los efectos de la indistinguibilidad de los fotones, es 
aquí donde ocurren los efectos cuánticos de la estadística de los fotones. Desde este punto de vista se puede afirmar 
que la onda electromagnética clásica es, en realidad, un efecto cuántico debido a la presencia de muchos fotones 
coherentes. 
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que cada modo nonnal del campo de radiación está cuantificado, de modo que su energía toma 
solamente los valores discretos 0 , hv, 2 hv, ..., ahora estamos considerando la radiación como un 
sistema de partículas idénticas (los cuantos del campo, es decir los fotones), cada una de las 
cuales tiene una energía hv. 

La emisión y absorción de fotones 

La naturaleza Bosónica de los fotones tiene importantes consecuencias sobre los fenómenos de 
interacción entre la radiación y la materia. Si bien para el estudio detallado de estos procesos 
hace falta la Electrodinámica Cuántica (EDC), que no vamos a tratar en estas notas, es posible 
deducir algunos aspectos fundamentales de la interacción radiación-materia en base a argumen¬ 
tos sencillos, sin necesidad de invocar la EDC. Eso es lo que vamos a hacer aquí. 

Emisión y absorción de fotones y la distribución de Planck 

Nuestra deducción de la estadística (15.64) del gas de fotones se basó en la distribución de Bose- 
Einstein, que se obtiene de la Mecánica Estadística considerando la distribución gran canónica 
de un sistema de Bosones en equilibrio a temperatura y volumen fijos, con el requerimiento adi¬ 
cional que el potencial químico del gas de fotones es nulo (porque la emisión y absorción de fo¬ 
tones por la materia determina el número total de fotones de la cavidad, de modo que éste no 
está determinado apriori). Vamos a mostrar ahora que se puede llegar a la (15.64) directamente, 
considerando explícitamente los procesos de emisión y absorción. 





-> 

hv=e b -e a 





-> 

hv=e b - e a 


e 


a 


w 


e 


a 


absorción emisión 

Fig. 15.8. Dos estados estacionarios particulares ip a y ip b de un átomo. En la transición 
a * b el átomo absorbe un fotón de energía e = hv = e b - e a . Viceversa, en la transición 
b —* a el átomo emite un fotón de la misma energía. 

Sean dos estados estacionarios particulares i¡> a y ip b de un átomo de las paredes de la cavidad, 
de energías 16 e a y e b , con e a < e b (Fig. 15.8). Entonces, en la transición a -» b el átomo ab¬ 
sorbe un fotón de energía 


e = hv = e b - e a 


(15.70) 


16 Aquí ay b indican el conjunto de los números cuánticos que identifican el estado. Para cada nivel de energía hay, 
en general, varios estados degenerados, pero nosotros estamos considerando ahora uno entre ellos en particular. 
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viceversa, en la transición b -* a el átomo emite un fotón de la misma energía. Ahora bien, hay 
muchos estados fotónicos (f)¡ con energía s (su número está dado por la (15.65)), que describen 
las diferentes direcciones de propagación y polarización del mismo. Nosotros estamos interesa¬ 
dos ahora en uno sólo de ellos, correspondiente a una particular dirección de propagación y una 
polarización determinada, que designaremos con (p k . Consideremos entonces una transición 
a -» b en la cual el átomo que inicialmente está en el estado ip a absorbe un fotón en el estado 
<j) k . La probabilidad P a -, b que en la unidad de tiempo ocurra esta transición está dada por 

Pa-*b = Pa P a,n k ->b,n k -\ (15.71) 

Aquí p a = e~P e “ /Z es la probabilidad que el átomo esté en el estado ip a , y 

T a,n k ^b,n k - 1 =l ) |2 (15.72) 


es la probabilidad de transición por unidad de tiempo desde el estado inicial i¡) a Q> n . al estado 

Y] % ^ 

final ipi/b,, _i del sistema (átomo + radiación ) y fi m k es el término del Hamiltoniano que des¬ 
cribe la interacción del átomo con el fotón (¡) k . Del mismo modo, la probabilidad P b ^ a que en la 
unidad de tiempo el átomo que se encuentra en el estado ip b emita un fotón en el estado (¡) k es 

P b^a = Pb P b,n k ^a,n k + 1 (15.73) 

donde p b = e~P gh / Z es la probabilidad que el átomo esté en el estado ip b , y 

T b,n k ^a,n k + 1 =l ) |2 (15.74) 


es la probabilidad de transición del estado inicial ip b <P n al estado final ip a <P n +1 . 

No conocemos H¡ nkk , pero claramente, es razonable suponer que debe ser de la forma 

W = Kt,k“k + h L,k“k (15.75) 


donde a k y á k son los operadores de creación y de aniquilación de un fotón en el estado (p k , cu¬ 
yos elementos de matriz (ecs. (14.32) y (14.35)) son, respectivamente 


a ky, 


V=(*n k -lA*n k ) = 


„ 4 - 72 , 1 . +1 


a 


t n k 
km,. 


-d’n,4l’nJ-p t +l 


(15.76) 


y h m k y su conjugado Hermitiano h¡ ntk son operadores 18 que actúan sobre las funciones de 
onda del átomo. La forma exacta de h int k y /z |nt k depende de las variables dinámicas que des¬ 
criben las cargas y corrientes eléctricas atómicas y de la función de onda cj) k del fotón que des¬ 
cribe los campos electromagnéticos asociados con el mismo; veremos que para nuestros fines no 
hace falta conocerla en detalle. Sustituyendo entonces (15.75) y (15.76) en (15.72) resulta 

rp a,n k ^*b,n k -1 — ^ ^ (p n k -l^k (p n k ) — n k \ b > hmt,k*P a ) (15.77) 


17 En la designación de <I> damos por sobreentendidos los números de ocupación que no cambian en la transición. 

IS Estos operadores describen la interacción de las cargas y corrientes atómicas con el campo electromagnético del 
fotón. 
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Haciendo las mismas sustituciones en (15.74) y observando que 

1 (Va’tíukQb) |2=l (Kuk^a^b) |2=l (Vb’hint.kVa) ^ (15.78) 

obtenemos 

T b,n k ^a,n k + 1 =l (Va’ÑnUkVb) ^ ( 0 n k+ l^l 0 n k ) |2 = ( n k + 1) ' 0/V W^a) |2 (15-79) 

Se puede notar que las (15.77) y (15.79) implican que las probabilidades de transición por uni¬ 
dad de tiempo de un proceso y su inverso son iguales, es decir 

T a,n k -*b,n k -1 = T b,n k -l^a,n k (15.80) 

lo cual significa que en los procesos de emisión y absorción se cumple el principio del balance 
detallado que mencionamos en el Capítulo 14. Esto es consecuencia de la Hermiticidad de h mtk , 
que se advierte en la (15.75). La presencia de los factores n k y n k +1 en (15.77) y (15.79) se 
debe a la naturaleza Bosónica del fotón. 

Ahora bien, si la radiación está en equilibrio térmico a la temperatura T se debe cumplir que 
= P b ^ a . Resulta por lo tanto que 

n k e^=n k +1 (15.81) 

Es importante observar que puesto que /),_/, y P b ^ a contienen ambos el factor (ip b J l \ nUk ip a ) 
éste se cancela al imponer la condición P a ^ b = P b ^ a > 1° cual implica que el mismo resultado 
(15.81) se obtiene al considerar las transiciones entre cualquier par de los estados atómicos de¬ 
generados de energías e a y e b que llevan a la emisión y absorción de un fotón (p k . 

De la (15.81) obtenemos de inmediato 

(,5 - 82) 

que es la (15.64). Vemos por lo tanto que la distribución de Planck es una consecuencia de tres 
hechos: (a) que la dinámica de los procesos de emisión y absorción de fotones cumple con el 
principio del balance detallado, (b) que los fotones son Bosones y (c) que hay equilibrio térmico. 
La física de la interacción entre átomos y fotones: emisión espontánea y estimulada 
Imaginemos un átomo que está efectuando una transición del estado ip a al estado ip b en la cual 
absorbe un fotón <¡> k de frecuencia v = (e b - e a )/h (o viceversa, una transición del estado i¡J b al 
estado ip a en la cual emite un fotón de esa frecuencia). En un instante dado durante la transición 
su estado está descripto por una función de onda de la forma 

W{t) = c a (t)ip a e~< llñ + c b (t)ip b e- ie » tlh (15.83) 

donde c a (t ) y c b (t) son ciertos coeficientes que dependen del tiempo (no nos interesa por el 
momento la normalización de W). La densidad de probabilidad es entonces proporcional a 

I W(t) l 2 =l c a (t) \ 2 \ip a \ 2 +\ c b (t) l 2 l ip b I 2 +2Re[cl(t)c b (t)ip* a ip b e- i{eh - ea)t/ñ ] (15.84) 
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y como se ve tiene una parte que oscila con la frecuencia v = (e b - e a )/h. Del mismo modo, si 
se calcula la corriente de probabilidad, se encuentra que contiene una parte que oscila con la 


misma frecuencia v. Puesto que la densidad de probabilidad y la corriente de probabilidad están 


relacionadas con las distribuciones espaciales de carga y corriente eléctrica del átomo, vemos 
que el estado W describe distribuciones de cargas y corrientes que oscilan con la misma frecuen¬ 
cia que lo hacen los campos electromagnéticos del fotón que está siendo absorbido (o emitido). 
Esta observación sugiere que el fenómeno de absorción (o de emisión) de radiación por parte de 


un átomo es análogo al aumento (o disminución) de la energía de un oscilador que está siendo 


forzado a oscilar por una fuerza aplicada, en resonancia con su frecuencia natural. Que el oscila¬ 
dor gane o pierda energía en este proceso, depende de la fase relativa entre sus oscilaciones y las 


de la fuerza excitadora. Esta es la imagen clásica de la absorción y emisión de radiación por 


parte de un átomo. De acuerdo con esta imagen es natural pensar que tanto la absorción como la 
emisión de radiación son procesos estimulados (o inducidos ) por la presencia del campo de ra¬ 
diación de la frecuencia que corresponde a la transición, y que su probabilidad debe ser propor¬ 
cional a la intensidad del mismo. 

Por otra parte los resultados anteriores (15.77) y (15.79) nos muestran que las probabilidades de 
transición por unidad de tiempo para la absorción y emisión de un fotón están dadas por 


n k I (ip b ,h mUk ip a ) I 2 absorción 

(n k + 1) I {ip b ,h ml k i¡) a ) I 2 emisión 



(15.85) 


Observemos que la intensidad del campo de radiación de frecuencia v (y dirección de propaga¬ 
ción y polarización correspondientes a (p k ) es proporcional a la energía del mismo, esto es, a 
n k hv. Vemos por lo tanto que la tasa de absorción es efectivamente proporcional a la intensidad 
del campo, tal como esperábamos. No ocurre lo mismo con la emisión, cuya tasa resulta propor¬ 
cional a n k + 1, lo que implica que el átomo puede emitir un fotón aún cuando la energía del 
campo de radiación de frecuencia v es nula. Este es un resultado puramente cuántico, que pro¬ 
viene del principio de incerteza, que prohíbe que el campo electromagnético sea estrictamente 
nulo 19 . Debido a ello, en ausencia de fotones el campo tiene igualmente fluctuaciones, y son es¬ 
tas fluctuaciones las que provocan la emisión del átomo cuando n k = 0. 

Es lógico entonces escribir T b n n +1 como la suma de dos contribuciones: la que se debe a la 
emisión espontánea provocada por las fluctuaciones de vacío del campo electromagnético y la 
emisión estimulada por los fotones (p k presentes, esto es 



emisión espontánea (15.85) 

emisión estimulada 


19 


Sin entrar en detalles, al tratar en forma cuántica el campo electromagnético, de acuerdo con los lineamientos 


expuestos en el Capítulo 14, se deben interpretar las ecuaciones de Maxwell como las ecuaciones de los operadores 
del campo. De resultas de ello las componentes de los campos E y B son operadores que satisfacen ciertas 
relaciones de conmutación, las cuales tienen como consecuencia (entre otras) que en el estado de vacío los campos 
eléctrico y magnético no son estrictamente nulos sino que fluctúan. 
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Por otra parte, como es obvio, la absorción es siempre estimulada. 

Tal como vimos anteriormente, el grueso del espectro de la radiación que proviene de fuentes 
naturales es incoherente pues corresponde a números de ocupación n k « 1. Por lo tanto la emi¬ 
sión estimulada es inobservable en esas condiciones: la emisión es en su totalidad espontánea. 
Veremos en breve bajo qué condiciones es posible observar la emisión estimulada y cuales son 
sus peculiares características. 

Los coeficientes de Einstein 

Las tasas según las cuales tienen lugar los procesos de absorción y emisión de radiación por 
parte de un átomo están determinadas, como acabamos de ver, por las probabilidades de transi¬ 
ción. Cabe observar, sin embargo, que las probabilidades de transición que consideramos en la 
Sección anterior se refieren a transiciones entre dos particulares estados atómicos ip a y ip b , ele¬ 
gidos entre los g a estados degenerados del nivel e a y los g b estados degenerados del nivel e b , 
asociadas con la emisión o absorción de un fotón de energía e = hv = e b - e a y con una determi¬ 
nada dirección de propagación y una dada polarización, correspondientes a un particular estado 
(p k . Por otra parte el estudioso está muchas veces interesado en conocer la intensidad de la radia¬ 
ción de una dada frecuencia v absorbida o emitida por un conjunto de muchos átomos debido a 
transiciones entre los niveles e a y e b , y no le importa saber en cuáles estados particulares ip a y 
ip b se encuentra cada átomo antes y después de la transición, ni tampoco en conocer la dirección 
de propagación y la polarización de los fotones involucrados. En este caso, las probabilidades de 
transición relevantes para la absorción (A) y la emisión (E) son 


(15.86) 



donde T e¡ es la probabilidad que el átomo que se encuentra en uno cualquiera de los estados 
del nivel e a efectúe una transición a uno cualquiera de los estados del nivel e b , sin que importe 
ni la dirección de propagación ni la polarización del fotón absorbido y es la probabilidad 

que el átomo que está en un estado del nivel e b efectúe una transición a un estado del nivel e a , 
cualesquiera sean la dirección de propagación y la polarización del fotón emitido. Estas proba¬ 
bilidades de transición se relacionan con las T an -*b, n -i y T b n ^ a n + i dadas por la (15.77) y la 
(15.79) por medio de 



(15.87) 


Aquí las sumas comprenden los g a estados degenerados del nivel e a y los g b estados degenera¬ 
dos del nivel e b , y la integral se extiende sobre todas las direcciones de propagación del fotón; el 


factor 2 toma en cuanta las dos polarizaciones del fotón, y la división por el factor de degenera¬ 
ción del estado inicial se debe a que hay que promediar sobre los diferentes estados iniciales de¬ 


generados puesto que todos ellos son igualmente probables. 

Es importante observar que debido a las transiciones originadas por los procesos radiativos y por 
las colisiones que ocurren en el medio, los estados atómicos no son estrictamente estacionarios, 
sino que tienen una “vida media” r¿, que depende del estado ip¿ y de la densidad y temperatura 
del medio, que determinan la frecuencia de las colisiones. Por este motivo las energías de los 
estados atómicos están indeterminadas dentro de una incerteza Ae¡ ~ hir ¡. Por lo tanto las líneas 
espectrales (sea de emisión como de absorción) no son estrictamente monocromáticas, pues los 


250 


15. Estadísticas Cuánticas 


fotones provenientes de distintos átomos tienen energías que difieren en cantidades del orden de 
Áe¡. Las integrales en (15.87) se tienen que extender entonces sobre un intervalo de energías del 
fotón del orden de Ae alrededor de e. Esto da lugar a un perfil de línea de ancho no nulo, des- 
cripto por una función f(v) cuya forma precisa depende los procesos que determinan la vida 
media del nivel. La función /(v) tiene un máximo en v = v 0 = £ / h y tiende rápidamente a cero 
al crecer I v - v 0 I; es habitual normalizarla de modo que 

00 

ff(v)dv = 1 (15.88) 

o 

Para tratar estos problemas es útil emplear los coeficientes de Einstein, que se relacionan con T A 
y T e , y que se pueden determinar a partir de magnitudes fácilmente medibles. Así, escribiremos 
que la probabilidad que en ausencia de radiación el átomo que está en un estado del nivel e b 
emita espontáneamente un fotón de frecuencia comprendida entre vy v + dv está dada por 

A b ^J(v)dv (15.89) 

donde A¡ ) _ a es el coeficiente de Einstein de emisión espontánea, y representa la probabilidad 
que en la unidad de tiempo el átomo efectúe espontáneamente la transición b -> a, cualquiera 
sea la frecuencia del fotón emitido. 

Supongamos ahora que el átomo está en presencia de radiación, cuya densidad espectral de 
energía es u(v). La presencia del campo de radiación cambia la situación anterior de dos mane¬ 
ras: (a) si el átomo está en el nivel e b la probabilidad de la transición b -> a cambia; (b) si el 
átomo está en el estado e a se torna posible la transición a -*> b por absorción de un fotón. Puesto 
que ambos efectos se deben al campo es lógico suponer que en el primer caso la probabilidad de 
emitir un fotón de frecuencia entre vy v + dv sea ahora 

[A b ^J{y) + B h ^ a u(v)f'(v)\dv (15.90) 

y en el segundo caso, la probabilidad de absorción sea 

B a ~, b u(v)f"(v)dv (15.91) 

Las nuevas funciones f'iy) y /"(v) tienen sus máximos en v = v 0 = e/h, tienden rápidamente 
a cero al crecer I v - v 0 I y están normalizadas del mismo modo que f(v). A priori, las tres fun¬ 
ciones/ /' y /" pueden ser diferentes, pero hay por lo menos un caso, aquél en que el átomo 
forma parte de un sistema en equilibrio termodinámico, en que / = /' = /" y por ese motivo de 
supone que son siempre idénticas. Las cantidades B b ^ a y B a ^ b son, respectivamente, los coe¬ 
ficientes de Einstein de emisión estimulada y de absorción. Los coeficientes de Einstein son 
constantes características de los niveles a y b, y se pueden calcular teóricamente en términos de 
los elementos de matriz de h intk usando las expresiones (15.87), o en su defecto se pueden de¬ 
terminar experimentalmente. 

Es fácil mostrar que entre los tres coeficientes de Einstein deben existir dos relaciones 20 , y por 
ese motivo una vez que se conoce uno de ellos se pueden calcular de inmediato los otros dos. En 


20 


Estas relaciones son consecuencia del principio del balance detallado y de la naturaleza Bosónica de los fotones. 
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efecto, consideremos la radiación en una cavidad y los átomos de sus paredes, en equilibrio ter- 
modinámico a la temperatura T. La relación entre el número de átomos en niveles ay b está dada 
por la distribución de Boltzmann 


^b_ _ 8b_ e -P(e h -e a ) _ Sb_ e -Pe 
8a 8a 


(15.92) 


El número de átomos que en la unidad de tiempo efectúan la transición b -» a (por emisión 
tanto espontánea como estimulada, y suponiendo / = /') resulta ser 

00 

N h j[A h ^ a + B b ^ a u(v)]f(v)dv = N b [A b ^ a + B b ^ a u(v 0 )] (15.93) 

o 

dado que m ( v ) varía muy lentamente con v. Del mismo modo, el número de átomos que en la 
unidad de tiempo efectúan la transición a -* b (suponiendo / = /" ) es 

00 

N a f B a ^ b u(v)f (v)dv = N a B a ^ b u(v 0 ) (15.94) 

o 

En el equilibrio termodinámico estos números deben ser iguales entre sí, de modo que 

i _ (jg 

N a A b ^ a + B b ^ a u(v 0 ) 


Usando la (15.92) y omitiendo el subíndice O, obtenemos 

8tL e -pe = B a ^ b u(v) 

8a + B b^a u ( v ) 


(15.96) 


La densidad espectral de energía radiante u(v,T ) está dada por la (15.66) y entonces tenemos 


J> b_ e ~Ps - 

§a 


B, 


Sithv 3 


’a—*b 


(eP £ -l )A b ^ a + B l 


8 Jthv 3 


(15.97) 


'b^>a 


Esta ecuación vale para toda T; en particular, para T -» oo se tiene que e^ E -> 1 y e > 1 y 
por lo tanto se debe tener que 

8b B b^a = 8a B a^b (15.98) 


Sustituyendo este resultado en la (15.97) se obtiene una identidad si 

8 Jthv 3 


A b -*a B b-*a ’ 


(15.99) 


Las relaciones (15.98) y (15.99) son las relaciones buscadas. Las obtuvimos para el caso del 
equilibrio termodinámico, pero puesto que se trata de relaciones entre constantes valen siempre. 
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Relación entre los coeficientes de Einstein y las magnitudes macroscópicas 

La variación de la intensidad de un haz de radiación de frecuencia v que atraviesa un espesor dx 
de un medio en el cual no tienen lugar procesos de emisión o de difusión es, por definición 

dl(v) = -Ku(v)dx (15.100) 

donde k es el coeficiente de absorción. La contribución del haz a la densidad de energía radiante 
está dada por 4 Jtl(v)/c, luego la variación de u(v) en el lapso dt = dxlc que emplea el haz en 
atravesar el espesor dx es 


du(v) = -Ku(v)cdt (15.100) 

Por otra parte, si la atenuación se debe a las transiciones a -» b de los átomos que se encuentran 
en el nivel e a , usando la definición (15.91) del coeficiente de Einstein B a ^ b , se debe tener que 

du(v) = -hvN a B a ^ b u(y)cdt (15.101) 

puesto que cada transición sustrae del haz un fotón de energía hv. Comparando (15.100) con 
(15.101) se obtiene la siguiente relación entre k y B a ^ b : 

K = —N a B a ^ b (15.102) 

c 


En nuestro argumento no tomamos en cuenta, sin embargo, que el haz de radiación estimula la 
emisión de los átomos del nivel e b , que se manifiesta como una absorción negativa. Luego para 
obtener el verdadero coeficiente de absorción tenemos que restar del miembro derecho de la 
(15.102) la cantidad hvN b B b ^ a / c. Por lo tanto el coeficiente de absorción correcto es 


hv 


hv 


K ' ~ (N a B a ^b ~ N b B h->a) ~ N a B a-*b\ 1 “ 


Nb B h-*a 

N a B a —* b , 


-fN a B a f 


1_ &l^l| (i5.io3) 

8b N a ) 


donde usamos la relación (15.98). Si los átomos están en equilibrio termodinámico a la tempe¬ 
ratura T obtenemos entonces 

k' = —N a B a ^ b ( \- e - hvlkT ) (15.104) 

c 

De momento que hemos incluido la emisión estimulada en el coeficiente de absorción, la emisi- 
vidad j está dada solamente por la emisión espontánea. Por lo tanto 

j = ^-N b A b ^ a hv (15.105) 

de donde resulta, recordando las relaciones entre los coeficientes de Einstein, la Ley de 
Kirchhoff-Planck: 


k' 


2hv 3 1 

c 2 ghvlkT _ j 


= B v (T) 


(15.106) 
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Si se ignora la emisión inducida, es decir, si se usa la expresión (15.102) del coeficiente de ab¬ 
sorción en lugar de la (15.104), se obtiene en lugar de la (15.106) la expresión 

7 O /i 

^ = —^-e~ hvlkT = B v (T) ( hv»kT ) (15.107) 

K C ¿ 

que es una aproximación de la Ley de Kirchhoff-Planck que vale para hv » kT , es decir en el 
límite en que los fotones se comportan como partículas clásicas. 

En el límite opuesto, cuando hv « kT , la emisión estimulada domina al punto que se puede ig¬ 
norar la emisión espontánea. En este caso, de la (15.104) resulta que 

i, (hv«kT) (15.108) 

ckT 

a partir de la cual se obtiene la bien conocida relación de Rayleigh-Jeans para la relación entre la 
emisividad y el coeficiente de absorción 

2 

^- = 2/m-j = B V (T) (hv«kT) (15.108) 

k' \ cJ 


en este límite de onda clásica la constante de Planck no interviene en B V (T). 

El efecto láser 

Consideremos la expresión (15.03) del coeficiente de absorción, que tiene validez general: 

(15.109) 

C \ Sb^a) 


En el equilibrio térmico tenemos que 

1 Nb _ ] _ g-hvlkT > q 

8b N a 


(15.110) 


de manera que k' > 0, lo que significa que la absorción domina sobre la emisión estimulada y 
por lo tanto el coeficiente de absorción es positivo, lo que significa que el haz de luz se atenúa al 
atravesar el medio. Sin embargo, en condiciones particulares, puede ocurrir que 

^^>1 (15.111) 

8b N a 


y en este caso se tiene que la emisión estimulada domina sobre la absorción: el coeficiente de 
absorción es entonces negativo. De resultas de ello el haz de luz se amplifica al atravesar el me¬ 
dio. Este fenómeno se denomina efecto láser 21 y toda fuente luminosa que aprovecha este efecto 
se denomina láser. Para que tenga lugar el efecto láser es preciso que las poblaciones de los ni¬ 
veles e b y e a que intervienen en la transición estén fuera del equilibrio térmico, de manera que 


21 El término láser es un acrónimo derivado de la frase “light amplification by stimulated emission of radiation” que 
describe en inglés el fenómeno de la amplificación de luz por emisión estimulada de radiación. 
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se pueda satisfacer la condición (15.110). Dicho en pocas palabras, la población del nivel de más 
alta energía ( e b ) debe ser significativamente mayor que la que correspondería al equilibrio, cir¬ 
cunstancia que en la jerga se denomina inversión de poblaciones. Hay diferentes técnicas 22 para 
conseguir este resultado, que no vamos a comentar porque el lector interesado las puede estudiar 
en la literatura especializada. Aquí nos limitaremos a describir brevemente las peculiares carac¬ 
terísticas de la luz láser y algunos interesantes y novedosos fenómenos a que da lugar, debido a 
los cuales el láser es un fértil campo de investigación sea pura como aplicada, además de tener 
numerosas aplicaciones (que es imposible tratar en el breve espacio de estas notas). 

El elemento esencial de todo láser es el medio activo, que puede ser sólido, líquido o gaseoso, 
donde se produce la inversión de poblaciones de los niveles de interés mediante una técnica 
oportuna (descargas eléctricas, bombeo óptico, reacciones químicas, etc.). El funcionamiento del 
láser puede ser tanto continuo (como el familiar láser de He-Ne) como pulsado. En ambos casos 
el medio activo tiene que estar colocado dentro de un resonador óptico (esencialmente un 
interferómetro de Fabry-Pérot) que selecciona una particular dirección, de modo que los fotones 
emitidos en esa dirección atraviesan el medio activo muchas veces estimulando la emisión de 
más fotones en cada pasada. De esa forma la población de un particular estado fotónico cp k se 
torna muy numerosa, y parte de la misma sale del dispositivo como un haz coherente (pues todos 
los fotones, que ocupan el mismo estado, están en fase entre si) y perfectamente colimado. Para 
muchas aplicaciones este haz es suficiente. Si se necesitan haces de mayor intensidad, se em¬ 
plean dispositivos más complejos que consisten de cadenas cuyo eslabón primario es el láser 
que acabamos de describir (que se llama oscilador ) el cual genera un haz que luego pasa sucesi¬ 
vamente por una serie de amplificadores, que son unidades cada una de las cuales contiene el 
medio activo, pero que a diferencia del oscilador no poseen resonador óptico. Como el nombre 
lo indica la función de estas unidades es amplificar el haz que sale del oscilador hasta alcanzar la 
intensidad que se desea. 

La radiación láser tiene características diferentes de la radiación que proviene de las fuentes na¬ 
turales (térmicas). Puesto que los fotones están en el mismo estado, son coherentes (mientras que 
los que provienen de fuentes térmicas son mayormente incoherentes). Esto hace que se com¬ 
porten como una onda electromagnética clásica (dado que n k » 1) lo que permite, entre otras 
cosas, observar con suma facilidad los fenómenos de interferencia y difracción, que normal¬ 
mente no se ponen en evidencia con la luz ordinaria a menos de tomar los particulares recaudos 
descriptos en los textos de Óptica a fin de evitar los efectos de la falta de coherencia de los foto¬ 
nes. La otra consecuencia de la coherencia es que los fotones de un rayo láser pueden interactuar 
colectivamente (y no sólo individualmente, como ocurre con la luz ordinaria) con la materia. 
Esto da lugar a fenómenos no lineales peculiares cuando la intensidad del haz láser es muy 
grande. 

Sin entrar en detalles, se puede decir que no hay límite (salvo en lo que se refiere al costo) a la 
potencia máxima del láser, y que enfocando el haz que produce es posible alcanzar enormes 
densidades de energía, imposibles de obtener con fuentes térmicas. Por ejemplo, con un láser 


22 Si bien la emisión estimulada fue reconocida por Einstein en su famoso trabajo de 1917, tuvieron que pasar cuatro 
décadas antes que Charles Townes, A. I. Schawlow, Nikolay Basov y Aleksandr Prokhorov demostraran la 
posibilidad de construir un láser óptico y que en 1960 Theodore Maiman desarrollara el primer láser de rubí. Desde 
entonces el progreso fue vertiginoso, y hoy existe una gran variedad de láser que se emplean en las más variadas 
aplicaciones, en las que se aprovecha la altísima colimación, la coherencia y la gran potencia de la radiación láser. 
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pulsado de vidrio-Nd, que emite luz de A = 1.06 pm, se han obtenido intensidades de 10 1N 
W/cm 2 en el foco 2 ’. Una sencilla cuenta muestra que esta intensidad corresponde a tener una 
densidad de energía de 200 cV/Á'' (1 Á 3 es el orden de magnitud del volumen de un átomo), que 
implica campos eléctricos del orden de 10 12 V/m. 

Sin necesidad de llegar a valores tan extremos, la interacción luz laser-materia con haces de gran 
intensidad se caracteriza porque se tiene siempre un gran número de fotones coherentes que inte¬ 
ractúan simultáneamente 24 sobre cada átomo iluminado por el haz. Ocurren entonces una varie¬ 
dad de fenómenos no lineales que no se observan con la luz ordinaria. Algunos de ellos han en¬ 
contrado importantes aplicaciones. Mencionaremos aquí el efecto fotoeléctrico multifotónico, la 
excitación multifotónica y la ionización multifotónica de niveles atómicos. 

En el efecto fotoeléctrico producido con luz ordinaria (incoherente) cada fotoelectrón absorbe un 
único fotón y se cumple la relación de Einstein (4.18) entre la frecuencia de la radiación y la 
energía cinética K max de los fotoelectrones 

K max= hv ~ w 0 (15.112) 

donde w 0 es la función trabajo, que es una característica del metal empleado como fotocátodo. 
Si v < v () = wq / h , los fotones no tienen energía suficiente para extraer fotoelectrones y no hay 
efecto fotoeléctrico, por grande que sea la intensidad (o sea el flujo de fotones) del haz de luz. 
Pero cuando los fotones son coherentes, y si la intensidad del haz es suficiente, se puede produ¬ 
cir el efecto fotoeléctrico multifotónico, en el cual cada fotoelectrón absorbe simultáneamente 
cierto número Wde fotones. En este caso, la (15.112) se debe modificar, y se tiene que 

K max = Ihí hv ~ w 0 (15.113) 

y entonces la frecuencia mínima necesaria no es v 0 sino v 05V - = v 0 / N = vv 0 1 Nh. Si el haz es 
muy intenso, A puede ser un número muy grande, y entonces el efecto fotoeléctrico ocurre para 
(casi) cualquier frecuencia, tal como sería de acuerdo con la teoría ondulatoria clásica del campo 
electromagnético. 

En la excitación multifotónica se provoca una transición atómica desde un nivel e a a un nivel e b 
de mayor energía usando fotones cuya energía es una fracción (la mitad, un tercio, ...) de la di¬ 
ferencia £ = e b - e a (Fig. 15.9a). En efecto, si W fotones idénticos (en el caso de la absorción N- 
fotónica) cada uno de los cuales tiene la energía hv^ = e/W llegan al átomo dentro de un inter¬ 
valo de tiempo suficientemente breve, el átomo responde como si fuesen un único fotón de la 
frecuencia correcta v = f//i = jVv jV y puede absorberlos, realizando la transición e a -» e b . Se 
trata de un proceso no lineal, y por lo tanto depende de la intensidad del haz láser. Si el nivel e b 
pertenece al continuo (Fig. 15.9b), el proceso se denomina ionización multifotónica. Además del 
interés puramente científico que presentan estos procesos, la excitación multifotónica tiene hoy 
una importante aplicación, la microscopía por fluorescencia multifotónica, que pennite obtener 
imágenes tridimensionales de muy alta resolución. 

23 Según la ley de Stefan-Boltzmann la radiancia de un cuerpo negro a la temperatura de la superficie del Sol (5700 
°K) es de 6xl0 3 W/cm 2 . Para alcanzar una intensidad de 10 lx W/cm 2 (pero distribuida en todo el espectro) con una 
fuente térmica, se precisa T > 2xl0 7 °K. Tan sólo una mínima fracción de la misma se encontraría en el intervalo de 
frecuencias correspondiente al entorno de A = 1.06 u m. 

24 Puesto que estos fotones son coherentes, se pueden describir por medio de un campo electromagnético clásico. 
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e b 


A 


-> 



->■ 

hVw= (e b ~ e a )/N 




A 




A 



->■ 

/iVjv= (e b - e a )HN 


e 


a 


(a) (b) 

Fig. 15.9. Procesos no lineales en la interacción entre la luz láser y los átomos: (a) excita¬ 
ción multifotónica; (b) ionización multifotónica. 


El gas perfecto de Fermiones 

De acuerdo con las ecs. (15.26) y (15.27), en un gas de Fermiones que no interactúan se tiene 

1 V. 1 


n- = 


,/?(£,■-fi) + 1 


= 2- 


,0(£,--íu) + \ 


(15.113) 


Introduciendo la densidad de estados por unidad de intervalo de energía f(e), la sumatoria de la 
segunda de las (15.113) se puede escribir como una integral 


= 2 


1 


1 / e P{e¡-H) + i 


i 


- (2i+1) l J^rr (2s+1)2,lV 


(2 ?n) 3/2 f e ll2 de 

Kh 2 ) ) + 1 

o 


(15.114) 


donde 5 es el spin, y usamos la expresión (9.19) de f(e), que se dedujo considerando los estados 
estacionarios de una partícula libre de masa m que se mueve en una caja de volumen V. El factor 
2s + 1 en (15.114) tiene en cuenta los diferentes estados de spin de las partículas. En lo sucesivo 
vamos a suponer que 5 = 1/2. Haciendo el cambio de variable z = fíe, la (15.114) se escribe 
como 


„, w (^) 

V h 2 ) 


3/2 


í 


z l,2 dz 

+ 1 


, wllüf) 

\ h 2 ) 


3/2 




í 


z ll2 dz 


e~ + e 




(15.115) 


Esta es entonces la ecuación que determina el potencial químico n, que como se ve, es función 
de T, Vy n. 

La energía de Fermi y las propiedades del gas de Fermiones en el cero absoluto 

Llamaremos energía de Fermi, por definición, a 

s F = fj,(T = Q,V,n) (15.116) 
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Vamos a ver ahora que e f debe ser positivo. En efecto, si se tuviera e f < 0, en el límite T —> 0 
tendríamos que p,/3 -> e F /3 -> y puesto que 


z ll2 dz 

o?' 4 . p! ! !' 



(15.117) 


se tendría que en ese límite el miembro derecho de la (15.115) se anularía, lo que es absurdo 
pues implicaría n = 0. Por lo tanto debe ser e F > 0. Tenemos entonces que para T = 0 los nú¬ 
meros de ocupación medios de los estados valen 


1 

'' ~ e P(£i-e F ) + \ 


rl si Ej < Ep 
|0 si £¡ > Ep 


(T = 0) 


(15.118) 


Este resultado tiene una clara interpretación: cuando T = 0 el gas de Fermiones está en el estado 
de mínima energía, en el cual las partículas ocupan los n estados ip¡ de menor energía, pues el 
principio de exclusión no permite que haya más de una partícula en cada estado (las diferentes 
orientaciones del spin están tenidas en cuenta por el factor 2s + 1 en la (15.114)). Estos estados 
son aquellos cuya energía es menor que e f , que es la energía del nivel más alto ocupado para 
T = 0. Esto permite determinar e f , pues para T = 0 la (15.114) (con s = 1/2) se reduce a 

(r-0) (15.H9) 

\h 2 ) { 3 \h 2 ) F 


de donde obtenemos 

2 m 18 Jt V ) 


(15.120) 


De la (15.120) se ve que la energía de Fermi es inversamente proporcional a la masa de las partí¬ 
culas, y proporcional a la potencia 3/2 de su concentración media n / V. 

La temperatura de Fermi se define por medio de e f = kT F y vale 


7V = 


h 2 


( 3 _ 1 \ 


2/3 


(15.121) 


y la velocidad de Fermi se define por medio de v F = ^2 e f /m . En la Tabla 15.1 damos los valo¬ 
res de e f , T f y v F para los electrones de conducción de algunos metales. 

Es fácil calcular el valor de la energía interna del gas de Fermiones en el cero absoluto: 

E = 2^ £ ¡ ñ i = AtcVÍ —^ j jV 3/2 d£ = —n£p = E 0 (T = 0) (15.122) 

i ' ^ ' o ^ 


Se puede mostrar (no lo haremos aquí) que la entropía del gas de Fermiones es nula en el cero 
absoluto, de acuerdo con la Tercera Ley de la Termodinámica. Usando la relación de Euler 25 


25 Ver el Capítulo 8 de Termodinámica e Introducción a la Mecánica Estadística. 
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E = TS-pV + un (15.123) 

es fácil calcular que en el cero absoluto la presión vale 

p = p F = ~£ F (T = 0) (15.124) 

Por lo tanto un gas de Fermiones tiene una presión no nula en el cero absoluto. 

Tabla 15.1: Valores de e F , T F y v F para los electrones de conducción de algunos metales. 


Elemento 

n/V (10 22 cm“ 3 ) 

e F (eV) 

7>(10 4 °K) 

v F (10 s cm/s) 

Li 

4.60 

4.7 

5.5 

1.30 

K 

1.34 

2.1 

2.4 

0.85 

Cu 

8.50 

7.0 

8.2 

1.56 

Au 

5.90 

5.5 

6.4 

1.39 


El potencial químico del gas de Fermiones 

Para calcular el potencial químico conviene expresar el número de partículas en función de la 
temperatura de Fenni. De la (15.121) obtenemos que 


8 Jt T7 /'2 mkT F \ 

= \ j 


3/2 


(15.125) 


Introduciendo esta expresión en la (15.115) y usando la (15.44) obtenemos entonces 





3/2 

Li 3 / 2 (-e^) 


(15.126) 


Invirtiendo esta relación se llega al resultado buscado, que se muestra en la Fig. 15.10, donde 
también representamos el valor clásico de p, que se obtiene en el límite « 1 y resulta ser 26 


^clásico 


--WlnlzV^ 7 ) 


3 / 2 ’ 




V Ú 2 ) 


(15.127) 


Para T «T F , la dependencia de u en T es muy débil. En esta región vale 27 la aproximación 


í 


z q ~ l dz 

e z -^ ± 1 


= +r{q )Li (+^)«^ 
q 


1 | 2 q(q - \)ji 2 \ 
3(3 ± 1)77 2 ) 


(■ r¡ ->oo) (15.128) 


Usando la (15.128) se obtiene fácilmente una expresión aproximada del potencial químico que 
vale para temperaturas bajas: 


26 La (15.127) difiere de la (15.51) debido al factor2 (que proviene del spin) en el argumento del logaritmo. 

27 El cálculo no es difícil pero algo engorroso y por eso no lo reproducimos aquí. El lector interesado lo puede 
encontrar en Landau y Lifschitz, Statistical Physics, Pergamon. 
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[i « e f 



(T «7» 


(15.129) 


Para T » T F el potencial químico tiende a Aí c iá s i c0 (ya para T > 2 T F la diferencia es muy pe¬ 
queña). 



T/T f 

Fig. 15.10. Potencial químico de un gas de Fermiones no interactuantes. Con la línea de 
trazos se muestra el valor clásico de ¡u, dado por la ec. (15.127). 

Los números medios de ocupación 

Es interesante calcular los números medios de ocupación, dados por la primera de las (15.113). 
La Fig. 15.11 muestra ñ(e) para diferentes temperaturas del gas. El gráfico es semilogarítmico 
pues así se aprecian mejor los apartamientos desde la distribución clásica de Maxwell- 
Boltzmann, que se representa por medio de rectas de pendiente 1 / kT. Se puede observar que a 
medida que se eleva la temperatura a partir de T = 0 °K, las partículas que estaban ocupando los 
estados de energía más alta, próxima a la energía de Fermi, se excitan y pasan desde los estados 
con s < ¡u a estados con s > p .. Para T « T F la mayor parte de estas excitaciones afectan partí¬ 
culas con energía cercana a s = p. 

Un gas de Fermiones se dice degenerado cuando su temperatura es inferior a la temperatura de 
Fermi; en el cero absoluto el gas está completamente degenerado. En un gas degenerado, el 
principio de exclusión, que no permite más de una partícula por estado, obliga las partículas a 
ocupar estados de mayor energía. Si la temperatura es muy alta (T » T h ), el sistema no es de¬ 
generado puesto que los números medios de ocupación de los estados son mucho menores que la 
unidad, y en ese caso el principio de exclusión tiene poco efecto en la práctica. El gas se com¬ 
porta entonces como un gas clásico. La transición entre el gas completamente degenerado en 
T = 0 y el gas clásico para T » T F es continua, pero la descripción matemática del gas en este 
rango de temperaturas es complicada. 
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e/e F 

Fig. 15.11. Número de ocupación medio en función de la energía del estado para un gas de 
Fermiones. Las curvas corresponden a T/T F = 0.01, 0.1, 0.5,1,1.5, 2 y 3. La línea hori¬ 
zontal de trazos corresponde a ñ(e = ¡u) = 1/2. Se puede observar que las curvas para 
T /T F = 1.5, 2 y 3 son casi rectas, lo que indica que el comportamiento del gas es clásico, 
como es de esperar porque ñ(e) «1. En cambio las curvas para T!T F = 0.01, 0.1 y 0.5 
difieren de las rectas clásicas donde ñ(e) ya no se puede despreciar frente a la unidad. 


Calor específico de un gas de Fermiones 

Para calcular el calor específico molar a volumen constante procederemos de manera semejante 
a como lo hicimos para el gas de Bosones. Calculamos primero la energía interna del gas 


E = 2 ^ e i n i 

i 


2 Y £ i s 3 - 3/2 f £ 3/2 d£ 

2j e P{Bi-n) + i 2 F J e^C-tO + l 

o 


(15.130) 


Usando la (15.44) obtenemos 


E = - 



3/2 

Li 5/2 (-^) 


(15.131) 


Usando ahora la definición (15.52) de c v , y la expresión 

/r 1 d[i _ 3 Li 3/2 (-e^ 1 ) 
kT~kdT~ 2 Li 1/2 (-e^) 


(15.132) 


que se obtiene derivando la (15.126) con respecto de T, obtenemos entonces 
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3/2 


2 LÍ 5/2( _í) 


3Li 3/2 (-^) 2 


2 Li 1/2 (-^) 


(15.133) 


En la Fig. 15.12 se puede apreciar el comportamiento del calor específico. 


c v /R 



T/T f 


Fig. 15.12. Calor específico molar a volumen constante de un gas de Fermiones. Se ve que 
c v es siempre menor que el valor clásico 3R/2 (representado por la línea de trazos) y para 
T <T F disminuye rápidamente con T. Para T « T F la dependencia de c v en T es lineal. 

Usando la (15.128) se obtienen fácilmente las expresiones aproximadas de E y c v para tempe¬ 
raturas bajas: 


£« 


3 

-ne F 


5ji 2 

1 +- 

12 



(T « T F ) 


(15.134) 


y 



(T « T F ) 


(15.135) 


El valor bajo de c v para T « T h se debe a que para un cambio de temperatura del orden de T, 
cabe esperar que la energía de una partícula aumente en kT. Pero para la mayoría de las partícu¬ 
las, con excepción de aquellas cuya energía es próxima a la energía de Fermi, los estados hacia 
los cuales se tendrían que excitar están ocupados (como se puede ver de la curva correspondiente 
a T = 0.01 T F en la Fig. 15.11). Por lo tanto, de acuerdo con el principio de exclusión, esas partí¬ 
culas no se pueden excitar. Sólo se pueden excitar las partículas que están en un intervalo de 
energía del orden de kT alrededor de e F , cuyo número es del orden de 

n exc ~f(s F )kT (15.136) 
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Por otra parte 


f(e F ) = 


AnV 
h 3 


(2 m) 3/2 f/ 2 


(15.137) 


Dividiendo esta ecuación por la (15.119) resulta entonces 

n 3 T 

'*exc _ _____ 

n 2 T f 


(15.138) 


de manera que la fracción de partículas que es excitada es del orden de T/T F . Puesto que cada 
una de ellas se lleva una energía del orden de kT, tenemos que 


E(T) ~E 0 +kT /z exc =~ n£ F 


-«a 


(T « T F ) 


(15.139) 


y el calor específico molar resulta entonces 


T 

c v « 3R — 
Tp 


( T«T F ) 


(15.140) 


Comparando la (15.139) y la (15.140) con los resultados exactos (15.134) y (15.135) vemos que 
nuestra estimación es correcta dentro de un factor del orden de la unidad. 

Discutiremos ahora brevemente un par de aplicaciones de la teoría que hemos desarrollado. 

El modelo de electrón libre de los metales 

En el estado sólido los átomos de un metal forman una red cristalina y los electrones de valencia, 
que están débilmente ligados, se pueden mover libremente en todo el cristal. Por eso cuando se 
aplica una diferencia de potencial eléctrico al metal, esos electrones producen una corriente, y 
por lo mismo se denominan electrones de conducción. Los electrones de conducción se mueven 
como partículas casi libres y su camino libre medio es muy grande (bajo condiciones adecuadas 
puede ser del orden del cm); en efecto, son poco perturbados por la presencia de los iones positi¬ 
vos y sólo sufren escasas colisiones con los demás electrones de conducción. 

Las razones de este comportamiento son de carácter cuántico, y las vamos a mencionar breve¬ 
mente sin dar demostraciones, que el lector puede encontrar si lo desea en la bibliografía. Si se 
resuelve la ecuación de Schródinger para un electrón en el potencial creado por una distribución 
periódica de iones positivos, se encuentra que la parte espacial de la función de onda de un 
electrón de energía e tiene la forma 

ipk ~ 8k( r ) e ‘ k r > k 2 =2m*£/ñ 2 (15.141) 


Aquí la función gfr) es \m factor de modulación que tiene la misma periodicidad espacial que 
la red, y el parámetro m = m "(k), que se denomina la masa efectiva del electrón, depende del 
vector de onda A: y de la geometría de la red. Si la red es infinita, resulta que no todos los valores 
de e están permitidos, sino que el espectro de energía consiste de una serie de tramos continuos, 
llamados bandas, separados casi siempre por intervalos finitos aunque puede ocurrir en determi¬ 
nados casos que dos bandas se superpongan parcialmente. Cada banda surge de la ruptura de la 
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degeneración de los infinitos niveles que consisten en que el electrón esté ligado en un particular 
nivel a un determinado ion de la red. Si el cristal no es infinito y consiste de N iones, entonces 
cada banda no es continua, sino que consiste de un conjunto de N niveles discretos, muy 
próximos entre sí. 

Lo anterior tiene una importante consecuencia en lo referente a la conductividad eléctrica del 
cristal. En efecto, para que la aplicación de un campo eléctrico externo de intensidad moderada 
produzca una corriente, es preciso que se pueda producir un movimiento neto de los electrones. 
Si todos los N niveles de una banda están ocupados, esto no es posible, pues los estados corres¬ 
pondientes a todos los valores de k están ocupados y no hay forma de que tenga lugar un flujo 
neto de electrones. Este es el caso de todas las bandas que se originan a partir de los niveles 
atómicos internos. La única banda que puede dar lugar a conducción de la electricidad es enton¬ 
ces la banda de valencia que surge a partir de los niveles atómicos más externos de los átomos, y 
las propiedades eléctricas del cristal dependen esencialmente de dos circunstancias: (1) cuántos 
electrones de valencia tiene cada átomo y (2) cuán separada en energía se encuentra la banda de 
valencia de la banda inmediatamente superior (originada a partir del primer nivel excitado 
atómico). Se pueden dar dos situaciones: (a) que la banda de valencia esté parcialmente llena, 
(b) que la banda de valencia esté completamente llena. En el primer caso, el cristal es un con¬ 
ductor, pues la aplicación de un campo eléctrico aún de muy baja intensidad puede excitar algu¬ 
nos electrones a estados vacíos de la banda, de modo tal que haya un flujo neto de electrones; 
esto es lo que ocurre con los elementos de valencia impar. En el segundo caso, que ocurre para 
los elementos de valencia par, si la separación en energía de la banda inmediatamente superior 
es considerable, el cristal es un aislante, pues la aplicación de un campo eléctrico moderado no 
produce corriente; esto sucede hasta que el campo es tan intenso que puede causar que algunos 
electrones pasen a la banda de energía más alta, que está vacía, en tal caso se produce la ruptura 
(o “breakdown”) del aislante. Pero si la banda inmediatamente superior está muy próxima en 
energía, si bien el cristal es aislante a T = 0, cuando la temperatura aumenta algunos electrones 
de la banda de valencia pueden pasar a la banda superior si kT « Ae, donde Ae es la separación 
en energía de las bandas; entonces el cristal puede conducir la electricidad. Este es el caso de los 
semiconductores. Finalmente, si la banda de valencia está llena pero hay una superposición par¬ 
cial con la banda inmediatamente superior, el cristal es conductor, este es el caso de algunos 
metales bivalentes. 

Volviendo ahora a los metales, un electrón de conducción en un cristal perfecto a 0 °K no está 
localizado, sino que se propaga libremente. Sólo las imperfecciones de la red (debidas a la pre¬ 
sencia de impurezas, los defectos de la red y las vibraciones de los iones) perturban la propaga¬ 
ción de los electrones. Además el fondo de carga positiva debido a los iones compensa en pro¬ 
medio las cargas de los electrones, y las interacciones residuales entre ellos tienen escasa im¬ 
portancia. Esto último es una consecuencia del principio de exclusión. En efecto, consideremos 
una colisión entre dos electrones que están inicialmente en los estados i¡J k y ip k *, y que después 
de la misma están en los estados xp k " y ip k ">. Esta colisión puede ocurrir solamente si los esta¬ 
dos ip k " y x¡) k "' no están ocupados, pues en el caso contrario este proceso está prohibido por el 
principio de Pauli. Por otra parte, la mayoría de los estados finales ip k „ y i¡> k que son energéti¬ 
camente posibles están ocupados. Por lo tanto estos procesos no ocurren y por este motivo es 
razonable pensar que en primera aproximación los electrones de conducción se comporten como 
si fuesen libres. 
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Esta es la justificación del modelo de electrón libre, que trata los electrones de conducción de un 
metal como un gas perfecto de Fermiones, ignorando su interacción con los iones positivos de la 
red, y las interacciones entre los electrones mismos. Pese a su simplicidad este modelo es muy 
útil pues explica en forma cualitativa muchas características de los metales. 

El modelo de electrón libre de los metales fue desarrollado originalmente por Paul Drude (1900) 
y Hendrik Antoon Lorentz (1905) usando la estadística clásica de Maxwell-Boltzmann. Se pudo 
así deducir la ley de Wiedemann-Franz (que da la dependencia con T del cociente entre las con¬ 
ductividades térmica y eléctrica), pero esto se debe a una coincidencia fortuita. Desde el inicio la 
teoría de Drude enfrentó un grave problema, pues tratados en forma clásica los electrones libres 
aportarían a la capacidad calorífica del metal una cantidad jk por electrón, que se tendría que 
sumar a la capacidad calorífica debida a las vibraciones de la red cristalina. De resultas de esto, a 
la temperatura ambiente se debería tener c v = (3 + 3/2)/? para un metal monovalente (con un 
solo electrón de conducción por átomo). Sin embargo los metales cumplen muy bien la Ley de 
Dulong y Petit, esto es c v =3R. 

La Mecánica Cuántica permite resolver esta dificultad. En efecto, lo correcto es tratar los elec¬ 
trones como Fermiones, lo cual fue hecho por primera vez por Sommerfeld en 1928. De la Tabla 
15.1 vemos que la temperatura de Fermi de los metales es de más de 10 4 °K. Por lo tanto, a tem¬ 
peraturas inferiores a la temperatura de fusión, se tiene que T « T F y entonces los electrones de 
conducción están degenerados. Podemos entonces usar la expresión aproximada (15.135) para 
calcular la contribución de los electrones a c v . Esta contribución es muy pequeña en compara¬ 
ción con la que proviene de las vibraciones de la red, lo cual explica porqué los metales cumplen 
muy bien la Ley de Dulong y Petit. 

La situación es diferente a temperatura muy baja. Cerca del cero absoluto, la capacidad calorí¬ 
fica de la red cristalina es proporcional a T 3 , mientras que la capacidad calorífica electrónica es 
proporcional a T. Por lo tanto a temperaturas suficientemente bajas domina la contribución elec¬ 
trónica, cosa que en efecto se verifica experimentalmente. 

El equilibrio de i as enanas blancas 

Se denomina enana blanca una estrella que al final de su evolución se ha transformado en un 
objeto muy denso, pequeño y poco luminoso. En una estrella normal (como el Sol) ocurren una 
serie de reacciones nucleares cuyo resultado neto es la fusión de 4 núcleos de hidrógeno para 
producir un núcleo de helio. En cada instante se producen gran número de tales reacciones, y en 
consecuencia se libera una enonne cantidad de energía. Esto es lo que mantiene a la estrella ca¬ 
liente e impide que se contraiga debido a la gravedad, manteniéndola en equilibrio con una den¬ 
sidad media baja (1.416 g cnL 3 para el Sol). Cuando la estrella agota su hidrógeno y pierde su 
fuente de energía, se enfría y finalmente se contrae. Cuando esto ocurre los átomos de helio se 
comprimen muy fuertemente. Los electrones que orbitan un átomo tienden a quedar más locali¬ 
zados y por lo tanto, debido al principio de incerteza, su energía cinética aumenta. Si la densidad 
aumenta suficientemente, la energía cinética de los electrones supera la energía potencial que 
tiende a ligarlos a los núcleos, y se convierten en partículas libres. De resultas de esto cada 
átomo de helio se transforma en una partícula a más dos electrones libres. Hay entonces una 
densidad enorme de electrones, que están relativamente fríos pues el combustible nuclear quedó 
agotado. Este es el estado de enana blanca. 

Para estimar la energía cinética por electrón cuando es una partícula libre podemos razonar del 
modo siguiente: cada electrón ocupa como mínimo un volumen del orden de 
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(15.142) 


donde Á B es la longitud de onda de Broglie. Si tenemos n electrones no podrán ocupar un volu¬ 
men menor que V = nvj .2 (dividimos por 2 pues dos electrones con spin opuesto pueden ocupar 
el mismo lugar) y entonces la energía cinética debe ser mayor que 


p 2 h 2 ( n ^ ?/3 
2 m e 2m e \2V) 


(15.143) 


Por otra parte la concentración de los electrones es ntV = Z/ d? donde d es la distancia entre los 
núcleos y Z su carga. Por lo tanto la energía cinética típica de un electrón es 


p 2 2 1/3 jr 2 Ú 2 Z 2/3 

2 m e m e d 2 


(15.144) 


el valor típico de la energía potencial es -2 Ze 2 / d . Por lo tanto el valor absoluto de la energía 
potencial es menor que la energía cinética si 
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(15.145) 


donde a 0 = ti 2 lm e e 2 es el radio de Bohr. Cuando esta desigualdad se satisface (a menos de un 
factor numérico del orden de la unidad) los electrones no pueden quedar ligados a los núcleos. A 
partir de la misma podemos estimar qué densidad debe tener la estrella para que los electrones 
queden libres. Sea m' la masa de la estrella por electrón. Puesto que al ionizarse cada átomo 
produce Z electrones, la masa de la estrella por núcleo es Zm !; por lo tanto la densidad tiene que 
cumplir (a menos de un factor numérico del orden de la unidad) la desigualdad 


P = 


Zm Z 2 
~d 3 “ 


m 


n 3 
a o 


(15.146) 


Para enanas blancas constituidas por helio Z = 2, luego m' = 2m H = 3.345 x 10 -24 g y por lo 
tanto la densidad tiene que ser mayor que 100 g/cm 1 . 

Las consideraciones anteriores presuponen que los electrones no son relativísticos, y para que 
esta sea una aproximación razonable se debe cumplir que v F « c, una restricción que implica 
que la densidad no debe superar unos 2xl0 6 g/cm 3 . Para estrellas cuya densidad se encuentra 
entre 100 y 2xl0 6 g/cm 3 , lo electrones se comportan como un gas degenerado de Fermi no rela- 
tivístico, en el cual las interacciones Coulombianas con los núcleos se pueden despreciar. 

La densidad de una típica enana blanca es de 10 6 g/cm 1 , la masa es de 10 J1 g, y la temperatura 
(interior) es de 10 7 °K. Para estos valores, T F « 2 x 10 9 °K, de modo que los electrones forman 
un gas de Fermi degenerado. La enana blanca más conocida es Sirio B, cuya masa de 2xl0 33 g 
es apenas menor que la del Sol; su radio es de 5600 km, algo más pequeño que el de la Tierra. 
Para ver cómo queda determinado el radio de una enana blanca, consideraremos primero el caso 
en que los electrones no son relativísticos. La presión en una estrella de masa M y radio R , de¬ 
bida a la gravitación, es del orden de 
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p G = - aGM 2 IR 4 


(15.147) 


donde G = 6.673 x 10 _8 dy cm 2 / g es la constante universal de la gravitación y a es un factor 
numérico del orden de la unidad. La estabilidad mecánica requiere que esta presión sea equili¬ 
brada por la presión del gas de electrones, dada por la ec. (15.124). El número de electrones es 
n = M/m' ( m' es la masa de la estrella por electrón) y V = 4 tcR 3 /3. Por lo tanto 


Pf = «' 


h 2 M 513 

m e m' 5,3 R 5 




(15.148) 


La condición de equilibrio p F + p G = 0 requiere entonces que 


M l/3 R = y 


ñ 2 

Gm e m' 513 



(15.149) 


Un cálculo preciso muestra que y «4.5. La (15.149) muestra que R ~ M~ 113 , y que es mucho 
menor que el radio de una estrella normal de la misma masa, de aquí la denominación de enana 
que se da a estas estrellas. 

La ley R ~ M~ l/3 se cumple mientras la densidad de equilibrio sea tal que los electrones no sean 
relativísticos. Pero p ~ MR~ 3 ~ M 2 , y por lo tanto si la masa es muy grande, al contraerse la es¬ 
trella los electrones se tornan relativísticos, lo cual invalida nuestro anterior cálculo del radio de 
equilibrio. Veamos qué sucede entonces, para lo cual conviene considerar el límite ultrarelati- 
vístico en el cual p 2 / m e » m e c 2 , y por lo tanto la energía cinética está dada por 

£ = ( p 2 c 2 + m 2 c 4 ) 1/2 - m e c 2 « cp (15.150) 


y no por la expresión clásica p 2 /2m e . En este caso no podemos usar la expresión (9.19) de la 
densidad de estados por unidad de intervalo de energía, y es preciso calcular nuevamente f(e) 
teniendo en cuanta que e = cp . Por medio de un razonamiento análogo al que nos permitió obte¬ 
ner la (9.19) es fácil verificar (lo dejamos como ejercicio para el lector) que ahora se obtiene la 
siguiente expresión de la densidad de estados por unidad de intervalo de energía: 


/(£) = 8 Jt 


V 


h 3 c 3 


(15.151) 


U sando la (15.151) resulta 


£p = he 


(jL±\ 
V 8jt v) 


1/3 


(15.152) 


La energía cinética total es entonces 



La presión se puede calcular a partir de la (15.123) y resulta 


(15.153) 
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(15.154) 


Para el caso de la estrella obtenemos 


p = y 


úc (M ) 473 

R 4 \m') 


y' -—( _ - _ ^ 4/3 

r 2 V 32 n 1 ) 


(15.155) 


Vemos entonces que cuando los electrones son ultrarelativísticos, la situación es completamente 
diferente al caso no relativístico, pues ahora la presión del gas de electrones es proporcional a 
R~ 4 , igual que la presión debida a la gravitación. La condición que ambas presiones sean iguales 
en valor absoluto detennina un valor crítico de la masa 


M 


■* 


( Shc \ 3/2 s _f 
\Gm ' 4/3 / ’ a 


(15.156) 


tal que si M < M' 1 ' el equilibrio es posible, pero si M > M* la presión del gas de electrones no 
puede nunca equilibrar la presión debida a la gravedad, y la contracción de la estrella continúa. 
Un cálculo exacto permite determinar que ó = 2.13. La masa crítica M" se denomina masa de 
Chandrasekhar 2S , y su valor es de aproximadamente 1.45 masas solares. Toda estrella cuya 
masa es menor que la masa de Chandrasekhar termina su existencia como una enana blanca 
estable, una vez que ha agotado su combustible nuclear 29 . 

Si M > M *, cuando la estrella ha agotado el hidrógeno colapsa bajo la acción de la gravedad. A 
medida que se produce dicho colapso, la densidad electrónica crece y por lo tanto el potencial 
químico de los electrones aumenta. Eventualmente, llega el momento en que se torna posible la 
reacción nuclear e~ + p -»• n + v e , en la cual un protón captura un electrón para transfonnarse en 
un neutrón y se emite un neutrino (esta reacción es la inversa del decaimiento ¡3 usual). De re¬ 
sultas de esto la estrella se convierte en una estrella de neutrones. La estabilidad de una estrella 
de neutrones 30 se puede estudiar del mismo modo que la de una enana blanca. Mientras la densi¬ 
dad no es demasiado elevada, se llega al equilibrio para un radio dado por una ley del tipo 
R ~ M _1/3 , como en el caso de una enana blanca. Pero si la masa es muy grande los neutrones 
se toman relativísticos y hay entonces un valor crítico de la masa, más allá del cual no hay equi¬ 
librio posible. Los cálculos se complican en este caso, pues es preciso usar la Teoría General de 
la Relatividad dado que el campo gravitatorio es muy intenso. No existe certidumbre sobre el 
valor crítico de la masa de una estrella de neutrones, pero se cree que es menor que 5, y más 
probablemente sea de alrededor de 3 masas solares. 

Las estrellas cuya masa es demasiado grande como para formar estrellas de neutrones estables 
no pueden detener el colapso y acaban su existencia como agujeros negros. 


28 Subrahmanyan Chandrasekhar, quien determinó dicho valor límite, obtuvo en 1983 el Premio Nobel por haber 
formulado la teoría de los estadios tardíos de la evolución de estrellas masivas. 

29 En particular, el Sol se convertirá en una enana blanca con un radio de unos 7000 km. 

10 En el curso del colapso la estrella puede explotar como una supernova, con lo cual pierde una fracción 
importante de su masa; en tal caso el residuo que se convierte en estrella de neutrones tiene una masa mucho menor 
de la que inicialmente tenía la estrella. 
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